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Abstract 

Heavy neutral gauge bosons, Z's, are predicted by many theoretical schemes of physics 
beyond the Standard Model, and intensive searches for their signatures will be performed 
at present and future high energy colliders. It is quite possible that Z's are heavy enough 
to lie beyond the discovery reach expected at the CERN Large Hadron Collider LHC, in 
which case only indirect signatures of Z' exchanges may occur at future colliders, through 
deviations of the measured cross sections from the Standard Model predictions. We here 
discuss in this context the foreseeable sensitivity to Z's of fermion-pair production cross 
sections at an e + e~ linear collider, especially as regards the potential of distinguishing 
different Z' models once such deviations are observed. Specifically, we assess the discovery 
and identification reaches on Z' gauge bosons pertinent to the E$, LR, ALR and SSM 
classes of models, that should be attained at the planned International Linear Collider 
(ILC). With the high experimental accuracies expected at the ILC, the discovery and the 
identification reaches on the Z' models under consideration could be increased substan- 
tially. In particular, the identification among the different models could be achieved for 
values of Z' masses in the discovery (but beyond the identification) reach of the LHC. 
An important role in enhancing such reaches is played by the electron (and possibly the 
positron) longitudinally polarized beams. Also, although the purely leptonic processes are 
experimentally cleaner, the measurements of c- and 6-quark pair production cross sections 
are found to carry important, and complementary, information on these searches. 
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1 Introduction 



Electroweak theories beyond the Standard Model (SM) based on spontaneously broken ex- 
tended gauge symmetries naturally envisage the existence of heavy, neutral, vector bosons 
Z 1 . The variety of the proposed Z' models is somewhat broad, and for definiteness in the 
sequel we shall focus on the so-called Z$ SM , Z'eri -^lr anc ^ -^alr models. Particular atten- 
tion has recently been devoted to the phenomenological properties and the search reaches 
on such scenarios, and in some sense we may consider these Z' models as representative 
of this New Physics (NP) sector [1,2]. 

A typical manifestation of the production of such states is represented by (narrow) 
peaks observed in the cross sections for processes among SM particles at high energy 
accelerators, for example in the invariant mass distributions for Drell-Yan dilepton pair 
production at the Fermilab Tevatron or at the CERN LHC hadronic colliders. Current 
experimental search limits on Mz> at 95% C.L., from Drell-Yan cross sections at the 
Tevatron, generally range in the interval 0.8-1 TeV, depending on the particular Z 1 model 
being tested [3]. Even higher 95% CL. limits, of the order of 1.14-1.4 TeV are obtained 
for the Z' x , Z^ R , and Z§ SM models, from electroweak high precision data [2]. 

Clearly, the eventual discovery of a peak should be supplemented by the verification 
of the spin-1 of the assumed underlying Z', vs. the alternative spin-2 and spin-0 hypothe- 
ses corresponding, e.g., to exchanges of a Randall-Sundrum graviton resonance [4] or a 
sneutrino [5]. This kind of analysis relies on appropriate angular differential distributions 
and/or angular asymmetries. Finally, once the spin-1 has been established, the particular 
Z' scenario pertinent to the observed signal should be identified, see, e.g., Refs. [6-15]. 
From studies of Drell-Yan processes at the LHC with a time-integrated luminosity of 
100 fb" 1 , it turns out that one can expect, at the 5-cr level, discovery limits on Mz> of 
the order of 4-4.5 TeV, spin-1 identification up to Mz' — 2.5-3 TeV and potential of 
distinction among the individual Z' models up to M z > — 2.1 TeV (95% C.L.). 

For masses above the direct search limits mentioned above, and LHC luminosity at 
the design value, access to Z' manifestations may be provided by indirect, virtual ex- 
change effects causing deviations of cross sections from the SM predictions, if Mz< is 
not excessively heavy. However, at the LHC, model identification from Drell-Yan dilep- 
ton mass distributions and forward-backward asymmetries may be problematic due to 
limited statistics [16]. 

An alternative resource for the observation of virtual heavy gauge boson exchanges 
should be represented by the next generation e + e~ International Linear Collider (ILC), 
with center of mass energy y^s = 0.5-1 TeV and typical time-integrated luminosities 
Ant ~ 0.5-1 ab _1 [17,18], and the really high precision measurements that will be possible 
there. Indeed, the baseline configuration envisages a very high electron beam polarization 
(larger than 80%). Also positron beam polarization, around 30%, might be initially 
obtainable and perhaps already available for physics. This polarization could be raised 
to about 60% or higher in the ultimate upgrade of the machine. The polarization option 
might represent an asset in order to enhance the discovery reaches and identification 
sensitivities on NP models of any kind [19], therefore also on Z' exchanges in interactions 
of SM particles. Previous analyses, based on various final state channels and possible 
experimental observables, show that sensitivities to quite high Z' masses could in principle 
be attained at the ILC (qualitatively, of the order of Mz' ~ (10 — 20) • ^/s for the highest 
planned luminosity, see, e.g., [1,20-24], and references therein). The ILC parameters have 
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recently been fixed in the Reference Design Report [17], so that it should be interesting 
to reconsider the identification of Z' models in the light of the numbers reported there. 

We will here focus on the fermion-antifermion production reactions at the polarized 
ILC: 

e + + e~^f + f, f = e,n,r,c,b. (1) 

As basic experimental observables for the Z' analysis, as an alternative to integrated 
observables like the total cross sections and/or angular- integrated asymmetries, we here 
choose the differential angular distributions for the above processes, that allow to exploit 
the information contained in the different portions of the final state phase space by a 
binned analysis. Particular emphasis will be given to the comparison between the cases 
of unpolarized and polarized initial beams, as regards the expected potential of ILC in 
identifying the Z' models of interest here, for Mz> values of the order of and beyond the 
limits accessible at the LHC. 

Indeed, concerning the Z' mass, there are two scenarios. The first one is represented 
by the interval in Mz> between the expected identification and discovery limits at the 
LHC: here, we can assume the Z ! to have already been discovered at some Mz> (but the 
model not identified), so that the model identification (or equivalently the determination 
of the coupling constants) could be performed at the ILC, based on the deviations of cross 
sections from the SM predictions for the determined Z' mass. For earlier attempts along 
this line see, e.g., Ref. [25]. The second mass range is above the LHC discovery limit and, 
here, with Mz> unknown, both discovery and identification reaches should be assessed for 
the ILC. 

In the following, in Sec. [2] we give a brief introduction to the different Z' models 
considered in the analysis, and give the corresponding leading order expressions of the 
polarized differential cross sections for processes (JTJ), mostly in order to establish the 
notations. In Sees. [3] and H] we present the results of our analysis for the discovery and 
identification reaches on the individual Z' models at the ILC; and finally, Sec. [5] contains 
some concluding remarks. 

2 Polarized observables and Z' models 

The analysis at the ILC is somewhat different from the corresponding studies of Drell-Yan 
processes at the LHC. Deviations of the various observables from SM predictions, such 
as cross sections and asymmetries, due to the interference of the SM amplitude with s- 
channel exchanges of the Z', graviton resonance G or sneutrino u, might be observed at the 
ILC. However, in the latter case, there is no interference of v with the SM exchanges [5]. 
Conversely, in the case of the spin-2 KK graviton exchange, the interference with the 
SM exchanges vanishes when one integrates over the full angular range [26], whereas for 
differential observables such interference survives. Nevertheless, it turns out [27] that the 
sensitivity at the ILC with \fs = 0.5 TeV and Ant = 500 ftr 1 to a KK graviton resonance 
in processes ([1]) is of the order 0.8 TeV (1.9 TeV) for the graviton coupling constant 
c = 0.01 (c = 0.1), which is well within the expectations for discovery and identification 
at the LHC. Accordingly, the KK excitation would have been either discovered or excluded 
by the time the ILC will be operating. Therefore, from the considerations above, one can 
conclude that for y/s < Mz>, as will be the case for the ILC, only the interference of the 
Z' amplitude with the SM one could be visible at the ILC in processes (pQ). Accordingly, it 
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might not be so indispensible to perform angular analyses such as those foreseen for Drell- 
Yan processes at the LHC, in order to differentiate the spin of the exchanged intermediate 
heavy quantum states, because only Z' should be able to lead to appreciable interference 
effects. 

The polarized differential cross section for the Bhabha process e + + e~ — > e + + e~, 
where 7 and Z can be exchanged also in the t-channel, can be written at leading order as 
(see, e.g., Refs. [28,29]): 

da(P-,P+) (1 + P~) (1 - P+) da R (1-P-)(1 + P+) da L 

dz 4 dz 4 dz 

(l + p-)(l + P+)dcr RL , f (1 - P-) (1 - P+) da LR , 

+ 4 dz + 4 dz ' [ ] 

with the decomposition 

dcr L _ da LL ^ da LR:S da R _ dcr RR ^ da RLjS 
dz dz dz ' dz dz dz 

In Eqs. (j2J) and the subscripts t and s denote helicity cross sections with SM 7 and 
Z exchanges in the corresponding channels, z = cos 8 and the subscripts L, R denote the 
respective helicities, P~ and P + denote the degrees of longitudinal polarization of the e~ 
and e + beams, respectively^ In terms of helicity amplitudes: 

dOLL _ 271T^, m , \ nee , r ee |2 d(T RR _ 27ra|^ , ,2 

dz " " |<*LL,, + <*LL,th d;z " fi |^RR, S + ^RR,*h 

dcrLR.t _ dgRLj _ 27ra^ m |rYee , 2 dcr LRiS _ da RLjS _ 2Tia 2 cm |rYee , 2 , ■ 
~dz~ ~ ~dz~ - I Gm -*I ' "^"^^"^^ |Glr ' s| • W 

According to the previous considerations the amplitudes with a,/? = L, R and 

i = s,t, are given by the sum of the SM 7, Z exchanges plus deviations representing the 
effect induced by a Z' boson: 



G?r„ = u - + -^br + , Gf 6 Tii = u - + t^tt + 



1 , G£) a (9 r 7 \ nee (1 , (9 R ) 2 (9*7 



'hh- « ' s s — M| s — M|, / ' t-My t-M\, 

1 Pr' 9l! \ „ee f 1 9 R 9l 9 R ' 9l 



Here, u,t = — s(l ± z)/2 (we are neglecting fermion masses), = — cot 2^ an d g R = 
tan 9y/ with 9y/ the electroweak mixing angle, whereas g' L and g R are characteristic of the 
particular Z' model. In the annihilation channel, below the Z' mass, the Z' interference 
with the SM will be destructive in the LL and RR cross sections, whereas it can be of 
either sign in the LR and RL cross sections. 

The polarized differential cross section for the leptonic channels e + e~ — > with 
I = H,t can be obtained directly from Eq. basically by dropping the t-channel con- 
tributions. The same is true, after some obvious substitutions, for the annihilations into 
cc and bb final states, in which case also the color (Nc) and QCD correction factors, 



^dn the recent review [19], the opposite sign convention for positron polarization was adopted. 
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C s ~ Nc [1 + a s /ir + 1.4 (a s /7r) 2 }, must be taken into account. The s-channel helicity 
amplitudes for the process with / ^ e,t can be written as [28]: 



n ef _ „ f QeQf 9 e a g f a 9a9'i \ r ef I Qe Qf 9a9p 9a 9 p 



Mf^ s-M*J' " aP ' s " \ s ^ s-M 2 z s-M% 

(6) 

where in the latter expression a ^ (3. 

As anticipated, the Z' models that will be considered in our analysis are the following: 

(i) The Z' scenarios originating from the exceptional group Eq spontaneous breaking 
are defined in terms of a mixing angle (3. The specific values (3 = 0, (3 = rr/2 and 
f3 = — arctan a/5/3, correspond to different E e breaking patterns and define the 
popular scenarios Z' x , Z'^ and Z' v , respectively. 

(ii) The left-right models, originating from the breaking of an SO (10) grand- unification 
symmetry, and where the corresponding Z£ R couples to a combination of right- 
handed and B — L neutral currents (B and L denote baryon and lepton currents), 
specified by a real parameter ctLR bounded by a/2/3 ^ a LR ^ a/2. The particular 
value «lr = v2 corresponds to a pure L-R symmetric model (LRS). 

(iii) The Z' ALR predicted by the 'alternative' left-right scenario. 

(iv) The so-called sequential Zg SM , where the couplings to fermions are the same as those 
of the SM Z. 

Detailed descriptions of these models, as well as the specific references, can be found, 
e. g., in Ref. [1]. All numerical values of the Z' couplings needed in Eq. ([5]) are collected, 
for example, in Table 1 of Ref. [14]. 



3 Discovery of Z' 



In the absence of available data, the assessment of the expected 'discovery reaches' on 
the various Z's needs the definition of a 'distance' between the NP model predictions and 
those of the SM for the basic observables that will be measured. The former predictions 
parametrically depend on the Z' mass and its corresponding coupling constants, while 
the latter ones are calculated using the parameters known from the SM fits. Such a 
comparison can be performed by a standard x 2 -like procedure. As anticipated in Sec. 1, 
we divide the full angular range into bins and identify the basic observables with the 
polarized differential angular distributions for processes (CQ), O = da(P~ , P + )/dz, in each 
bin. Correspondingly, the relevant x 2 can symbolically be defined as: 



/ {P-, P+} bins (OUbmJ 



Notice that not only the different beam longitudinal polarizations, but eventually 
also the various processes '/' in Eq. ([1]) are combined in the definition ([7]). Here, one 
assumes to have produced a set of 'data', for example by using the dynamics specified by 
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Figure 1: Discovery reaches on Z' models obtained from combined analysis of the unpo- 
larized and polarized processes (pQ) (95% C.L.) at the ILC with ^/s = 0.5 TeV (1 TeV) and 
Ant — 500 fb _1 (1000 fb" 1 ), compared to the results expected from Drell-Yan processes 
at the LHC at the 5-er level [14]. Three options of polarization are considered at the ILC: 
unpolarized beams, P~ = P + = 0; polarized electron beam, |P~| = 0.8; both beams 
polarized, |P~| = 0.8 and \P + \ = 0.6. 
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a given Z' model, and SO in the denominator denotes the corresponding 'experimental' 
uncertainty on O, combining statistical and, if possible, systematical ones. According to 
the previous considerations, the x 2 ; besides the number of degrees of freedom, is basically 
a function of the chosen Z' model parameters. In particular, if the coupling constants are 
fixed at specific values, it will depend solely on the Z' mass, and we vary this parameter. 
The discovery sensitivity to the Z' under consideration can in this case be identified as 
the limiting value of Mz> for which the value of x 2 {Mz>) has the probability needed for 
exclusion of the SM at a desired confidence level (in what follows, we shall impose 95% 
C.L.). In the cases where cos/3- or aLR-dependent couplings are considered, SM exclusion 
regions can be defined analogously. 

To derive the expected 'discovery' limits on Z' models at the ILC, for the 'annihi- 
lation' channels in Eq. (JT]), with f ^ e,t, we restrict ourselves to combining in Eq. (jTJ) 
the (P~,P + ) = (\P~\,— \P + \) and (— |P~|,|P + |) beam polarization configurations, that 
are the predominant ones. For the Bhabha process, / = e, we combine in (J7]) the cross 
sections with all four possible polarization configurations, i.e., (P~,P + ) = (|P~|, — |P + |), 
(-|P-|, |P+|), (|P-|, |P+|), (-|P-|,-|P+|). Numerically, following the ILC Design Re- 
port [17], we take for the electron beam |P _ | = 0.8. For the positron beam, |P + | = 0.3 is 
discussed as possibly available 'free of charge' already in the ILC initial running conditions. 
However, such a small positron polarization will turn out not to affect our evaluated dis- 
covery and identification reaches on Z's considerably. We shall therefore present numerical 
results for two cases, unpolarized positron beam |P + | = 0, and |P + | = 0.6 representing 
the 'ultimate' upgrade. 

Regarding the ILC energy and the time-integrated luminosity (which, for simplicity, we 
assume to be equally distributed among the different polarization configurations defined 
above), still according to Ref. [17], we will give explicit numerical results for cm. energy 
y/s = 0.5 TeV with time-integrated luminosity C m t = 500 fb~ , and for the 'ultimate' 
upgrade values yfs = 1.0 TeV with £ int = 1000 fb -1 . The assumed final state identification 
efficiencies governing, together with the luminosity, the expected statistical uncertainties, 
are: 100% for e + e~ pairs; 95% for l + l~ events (I = /i,t); 35% and 60% for cc and bb, 
respectively [17,18]. 

As for the major systematic uncertainties, they originate from errors on beam polar- 
izations, on the time-integrated luminosity, and the final-state reconstruction and energy 
efficiencies. For the longitudinal polarizations, we adopt the values <5P _ /P~ = 5P + /P + = 
0.25%, rather ambitious, especially as far as P + is concerned, but strictly needed for con- 
ducting the planned measurements at the permille level, see, e.g., Refs. [30-32] H As 
regards the other systematic uncertainties mentioned above, we assume for the combina- 
tion the (perhaps conservative) lumpsum value of 0.5%. The systematic uncertainties are 
included using the covariance matrix approach [33-35]. 

Concerning the theoretical inputs, for the SM amplitudes we use the effective Born 
approximation [36] vertices, with m top = 175 GeV and mn = 120 GeV. The numerically 
dominant 0(cx) QED corrections are generated by initial-state radiation, for both Bhabha 
scattering and the annihilation processes in ([1]). They are accounted for by a structure 
function approach including both hard and soft photon emission [37], and by a flux factor 
method [38], respectively. Effects of radiative flux return to the s-channel Z exchange are 

2 For simplicity we here assume equal precisions on beam polarizations. Clearly, much larger systematic 
errors on P + could partially spoil the advantages expected from the availability of also positron beam 
polarization. 
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Figure 2: Relative deviation of the unpolarized differential cross sections from the ALR- 
model prediction A for the process e + e~ — ► bb as a function of cos 9 for the SM and 
considered Z' models with Mz> = 2.5 TeV at the ILC with y/s = 0.5 TeV and C- mt = 
500 fb _1 . The error bars are statistical uncertainties at the 1-a level. 

minimized by the cut A = Er 1 /E\ )Bam < 1 — M§/s on the radiated photon energy, with 
A = 0.9. In this way, only interactions that occur close to the nominal collider energy 
are included in the analysis and, accordingly, the sensitivity to the manifestations of the 
searched-for nonstandard physics can be optimized. By numerical studies based on the 
ZFITTER code [39], other QED effects such as final-state and initial-final state emission 
are found, in the processes e + e~ — > (I = fi,r) and e + e~ — > qq (q = c,b), to be 
numerically unimportant for the chosen kinematical cuts. Finally, correlations between 
the different polarized cross sections (but not between the individual angular bins) are 
taken into account in the derivation of the numerical results, that we present in Fig. [TJ 
The figure includes a comparison with the discovery potential of the LHC with luminosity 
100 fb" 1 , from the Drell-Yan processes pp — > + X (I — e, fj) (at the 5-cr level). These 
values provide a representative overview of the sensitivities of the reach in Mz> on the 
planned energy and luminosity, as well as on beam polarization. 

4 Distinction of Z' models 

Basically, in the previous subsection we have assessed the extent to which Z' models can 
give values of e + e" differential cross sections that can exclude the SM hypothesis to a 
prescribed CL. Such 'discovery reaches' are represented by upper limits on Z' masses, for 
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which the observable deviations between the corresponding Z' models and SM predictions 
are sufficiently large compared to the foreseeable experimental uncertainties on the cross 
sections at the ILC. 

However, since different models can give rise to similar deviations, we would like to 
determine the ILC potential of identifying, among the various competing possibilities, the 
source of a deviation, should it be effectively observed. These ID-limits should obviously 
be expected to lie below the corresponding ILC discovery reaches and, for an approximate 
but relatively simple assessment, we adapt the naive x 2 -like procedure applied in the 
previous subsection. 

To this purpose, we start by defining a 'distance' between pairs of Z 1 models, i and 
j with i,j denoting any of the SSM, SM, ALR, LRS, ip, r], x, but i ^ j. We assume 
for example model i to be the 'true' model, namely, we consider 'data' sets obtained 
from the dynamics i, with corresponding 'experimental' uncertainties, compatible with 
the expected 'true' experimental data. The assessment of its distinguishability from a j 
model, that we call 'tested' model, can be performed by a x 2 comparison analogous to 
(EE]), with the x 2 defined as: 

x 2 (o)„, = j2 e e [g(z ;!;? L (s) 

/ {P~, P+} bins l<WbinJ 

As an illustration, the angular behavior of the deviations in the numerator of Eq. (jSj) 
for the unpolarized annihilation e + e~ — > bb is depicted in Fig. [21 for the case where 
the 'true' model is % — ALR, with M z > = 2.5 TeV for all models, at the ILC with 
y/s = 0.5 TeV and C- mt = 500 fb" 1 (actually, in this figure, A is the relative deviation, 
A = da(Z' ALR )/da(Z>) - 1). 

Basically, considering that the ILC will start when the LHC will already be operating 
at the design energy and luminosity, as anticipated previously, we can envisage two cases 
requiring somewhat different strategies. 

4.1 Z' mass known 

In the first case we assume that the Z' mass is already measured at the LHC, but perhaps 
not 'identified' there, and the value is within the ILC discovery reaches for both models 
i and j. In this case one should set M z > = M z >. = Mz> in Eq. ([8]) and, accordingly, the 

1 3 

X 2 becomes a function of only the Z[ and Z'^ coupling constants. If both the Z[ and 
couplings are fixed numerically, like in the example of Fig. [2j distinguishability can be 
assessed by varying M Z ', up to the point where the x% reaches the critical value suitable 
for exclusion of the 'tested' model j by the 'true' model i at the desired confidence level. 

If the above mentioned couplings are, instead, the (3- or «lr- dependent ones, 'confu- 
sion' domains between 'true' and 'tested' models can analogously be determined by means 
of Eq. ([8]) in the model parameter plane (cos (3, «lr) for fixed values of M Z ' = M z < = Mz>- 

1 J 

By definition, in these 'confusion' domains, that depend on the actual value assumed for 
M Z ', the cross sections corresponding to definite values of (3 and cklr cannot be distin- 
guished from each other at the desired confidence level. Correspondingly, the 'complemen- 
tary' regions in the above mentioned parameter plane can define the 'resolution' domain 
of the 'tested' model by the the 'true' model hypothesis, and determine in this way the 
identification limit on the latter. 
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e + e-^e + e~, M v = 4.5 TeV e + e~^l + r, M v = 4.5 TeV 




cos/? cos/? 



Figure 3: Regions of resolution (red hatched) of Eq vs. left-right LR models in the model 
parameter plane (cos /3, «lr) for M z > = 4.5 TeV obtained with Eq. (JSJ) from the processes 
with different final states e + e~ — > e + e~, (I = /x, r) cc, 66, at yfs = 0.5 TeV, C- in t = 
500 fb _1 . The role of polarization is demonstrated. Using also a polarized electron (and 
positron) beam, the resolution region would enlarge to include also the blue hatched 
(plus the shaded, unhatched) regions, but not the white ones. Those would remain as 
'confusion' regions. 
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Figure 4: Similar to Fig. [3] but for Mz> = 3.5 TeV. 



As an illustrative example of application of the ID-criteria exposed above, we evaluate 
the 'resolution' regions between Eq ('true') and left-right LR ('tested') models in the 
plane (cos/3,aLR) for different values of M z >. Figures [SJEl show the regions of 'resolution' 
obtained from the processes e + e~ — > e + e~, (I = fi, r), cc and bb, for Mz> =4.5 TeV, 
3.5 TeV and 2.5 TeV at y/s = 0.5 TeV and C int = 500 fb _1 , and for different values 
of beam polarization. Notice that, in these figures, the horizontal axis includes also the 
values of (3 specific of the x, i> and rj models, while the vertical axis includes the value of 
«lr representative of the LRS mo del@ 

Figures [3HS] clearly demonstrate the complementary roles of the processes with different 
final states, in particular, that the process e + e~ bb can potentially be the most efficient 
one in distinguishing Eq and left-right models from each other (it provides the largest 
resolution domains). Conversely, the purely leptonic processes, / = e, fi, r in (0Q), turn out 
to determine much less extended 'resolution' areas, in particular they cannot discriminate 



3 Actually, we should recall that the point (cos/3, q;lr) = (1, corresponds to the \ model which, 

indeed, exists in both classes of Z 1 models. Therefore, since one would be testing the \ model against 
itself, that point must not be considered in an analysis based on Eq. ©. 
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Figure 5: Similar to Fig. [3] but for M z > = 2.5 TeV. 

the LR model from the E e models at (cos/3, «lr) = (±1, v^/3) an d (1/4) a/3/2) (for any 
M z , > v^). 

Also, as can be seen from these figures, the leptonic processes are found to provide 
'confusion' domains (white) located in the 'central' part of the plane (cos/3,«lr), around 
(1/4, y/3/2), whereas the processes into qq final states exhibit the opposite feature. There- 
fore, as shown in Fig. [61 the combination of all processes / is expected to dramatically 
reduce the 'confusion' area in the above mentioned plane and to determine the largest 
possible domain in which the considered Z' models can be mutually distinguished from 
one another. The substantial role of electron polarization and, to a somewhat lesser ex- 
tent, of positron polarization in shrinking 'confusion' domains, leading to enlarged model 
'resolution' domains, can also be seen in these panels. Combining all processes, and with 
both beams polarized, a 'confusion' turns out to persist only in the minute corners shown 
in Fig. [HI for Mz> = 4.5 GeV, and nothing at the lower masses. 

It is interesting to compare these resolution regions with the corresponding ones re- 
sulting from the assumed Z' discovery in the Drell-Yan process at the LHC, shown in 
the lower-right panel of Fig. [6j This figure shows that, at the LHC, for the discovery of 
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>ff,M z > = 4.5 TeV 



pp->l + r+X; LHC@100fb _1 

5 




Figure 6: Left: Resolution domains at M z > = 4.5, exploiting all final states. Right: 
Corresponding resolution domains for Drell-Yan production at the LHC, for the three 
values of Mz>- 



a 4.5 TeV Z' the corresponding resolution region is found to cover only a narrow strip, 
1.4 ^ «lr and —0.5 ^ cos/? ^ 0.8. Even for M z > = 2.5 TeV, at the LHC, there are only 
modest hyperbola-like strips at | cos/5| ^ 0.5 where models can be distinguished. 

We now continue the above analysis in a somewhat different direction, namely, we 
wish to determine the limiting value of Mz> up to which a particular cos ^-dependent E 6 
model, assumed to be true, can be identified at the ILC in the sense that all the other, 
potentially competing, Z' models can be excluded. The results from this analysis are 
shown in Figs. EHHJ 

Figure [7J exhibits the exclusion limits vs. (3 on the models LR, ALR, SSM and SM 
(recall that exclusion of the SM determines the discovery reaches), once a Eq model is 
'true' (all processes combined). For the LR 'tested' model, the corresponding curve in 
Fig. [7J is obtained, for each (3, by varying «lr m the full allowed range, which gives LR 
exclusion limits Mz^a^n), and choosing the minimum value of such Mz'S (in this way the 
whole class of LR models, as well as the LRS, are excluded). The solid line labelled 'SM' 
represents the discovery reach, i.e., the Z' mass up to which the SM can be excluded. 
The overall identification range is shown as the shaded (yellow) region. One can see 
that, in this case, the identification of the class of E$ models considered here is basically 
determined by the exclusion of the class of the LR models and, for the 'central' values 
of cos/?, by the SM (i.e., by the discovery reach). Consequently the ID-limit is, in a 
(somewhat broad) range around cos/3 = 0, essentially identical to the discovery limit, 
whereas it is substantially smaller in the two intervals close to | cos /3| = 1. Figure [7J 
shows that, numerically, for | cos(3\ < 0.9 the ID-limit is as large as M%} ~ 3 - 4 TeV, 
and for cos f3 near ±1 E e models become more and more difficult to distinguish from the 
competitor ones. The right panel of Fig. [TJ shows the corresponding identification reaches 
for the polarized case, \P~\ = 0.8 and |P + | = 0.6, and the quantitative improvements 
that can be achieved in this case. 

Similarly, the identification limits on LR models vs. the parameter cklr can be read 
off from Fig. [SJ The curve labelled as 'E 6 ' is obtained by a procedure analogous to the 
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Figure 7: Left panel: E6 identification reaches on Mz> at 95% CL. obtained from 
combination of all unpolarized processes e + e _ — > f f at ^=0.5 TeV and £i nt =500 fb -1 . 
The Eq model is assumed to be 'true' while the others (SSM, LR, ALR, SSM and SM) 
are taken as tested models. The identification range is indicated as the shaded (yellow) 
area. Right panel: Similar, but for the polarized processes. 
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Figure 8: Similar to Fig. for LR models. 



curve 'LR' in Fig. [3, and the solid curve 'SM' represents the exclusion limits of the SM 
(hence the discovery reaches). In this case, the identification of the class of LR Z's turns 
out to be determined basically by the exclusion of the class of E e models, generally not 
so much below the discovery limit for all values of olr. On the other hand, the figure 
shows rather high identification limits, of the order of M|? ~ 3.0 — 4.6 TeV in the range, 
say, 0.9 ^ oxr ^ y/2, whereas they substantially decrease for smaller olr. 

We can conclude, from Figs.[7]and[8], that the identification reach at the ILC, already at 
\fs = 0.5 TeV and £ int =500 fb -1 , exceeds the corresponding discovery reach at the LHC. 
In fact, the full integrated luminosity considered here might be not quite indispensible for 
this identification. In Tabled] we show the required integrated luminosity, at the two ILC 
energies of 0.5 and 1 TeV, for the identification of these different models, realized as a Z' 
at 2.5, 3.5 or 4.5 TeV (within the discovery reach of the LHC). 

Finally, in Fig. [9] we summarize the information, of a similar kind as represented in 
Figs. [7] and [HI relevant to the cases where the ALR model or the SSM model is assumed 
'true' (upper and lower panels, respectively). As usual, the figure shows the limiting value 
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Figure 9: Exclusion reaches on Mz> at 95% CL. obtained from combination of all processes 
e + e' -> // at ^=0.5 (1.0) TeV and £ int =500 (1000) fb" 1 in the case when the ALR 
model (top panel) or SSM model (bottom) is assumed to be the 'true' models while the 
others are taken as tested models. The unpolarized and polarized cases are compared. 



Table 1: Required integrated luminosity, £ int [fb _1 ], at the two energies yfs = 0.5 and 
1 TeV and with polarized beams, required for model identification. Three mass values, 
Mz> = 2.5 TeV, 3.5 TeV and 4.5 TeV, assumed determined at the LHC, are considered. 





E6 

2.5 3.5 4.5 


LR 

2.5 3.5 4.5 


ALR 
2.5 3.5 4.5 


SSM 
2.5 3.5 4.5 


0.5 TeV 
1 TeV 


49.6 225 785 
9.1 41.6 125 


51.0 241 944 
9.4 42.3 128 


2.3 9.4 27.1 
0.4 2.0 5.9 


2.1 8.5 24.0 
0.4 1.8 5.3 



of Mz> at which the other Z' models can be excluded, the SM being one of them. The 
two energies 0.5 TeV and 1 TeV are considered, the different processes ([T]) are combined 
in the x 2 , and the unpolarized and polarized cases (\P~ \ = 0.8, \P + \ = 0.6) are compared. 
The entries for the LR and E$ models in this figure refer to the worst case, i.e., similar 
to the procedures adopted for the 'LR' and 'E$ curves in Figs. [7] and El adopting the 
lowest value of Mz> as (3 and q;lr are varied, in order to represent a whole class of models. 
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Specifically, in Fig. [9], one can easily read off the identification reaches for ALR and SSM 
models at 0.5 TeV: Ml 1 ? and Ml 1 ? are both 8-9 TeV if polarization is available. In a 

^ALR ^SSM 

sense, the ALR and SSM models are the most 'orthogonal' ones since, if either of them is 
assumed 'true', the other one can be excluded up to a really high value of M Z i. 

4.2 Z' mass not known 

The second kind of situation is met in the case where the Z' mass cannot be known a priori, 
e.g., the Z' is too heavy to be discovered at the LHC [say, M z > >4-5 TeV], but deviations 
from the SM predictions can still be observed at the ILC Actually, models with different 
Z' masses and coupling constants can in principle be the source of a deviation from the 
SM predictions observed at the ILC. With the coupling constants held fixed numerically 
at the theoretical values pertinent to the Z[ and Z', models under consideration, the xjj 
of Eq. (jHJ) becomes a function of the two masses, M z > and M z >., both assumed to lie 
in the respective ILC discovery ranges. In this case, one can derive a contour in the 
two-dimensional (Mz>, M z i) plane where to each value of Mz< is associated a value Mz 1 

such that, for all M z >. > M z >_, the value of x% i* 1 © is consistent with 'confusion' of i 
and j at the desired confidence level. The region encircled by such a contour will be the 
'confusion' (or 'no distinction') domain between the 'true' model i and the 'tested' model 
j and correspondingly, in the complementary domain the hypothesis j could be excluded 
if % is assumed to be 'true'. We refer to this latter, complementary, region as 'resolution' 
region. 

One can iterate this procedure and generate pairwise 'confusion' and 'exclusion' regions 
in the two-dimensional planes of parameters for all models j ^ i. As will be illustrated 
graphically in the remaining part of the paper, a common feature of such 'exclusion' 
regions is that the relevant contours admit, for each j (and obviously fixed i) a minimum 
value M^) such that, for any value Mz'. < M%) , the 'tested' model j can be excluded 
regardless of Mz 1 - We finally assume, as identification limit on the i model at ILC, the 

smallest of the values M~?for j ^ i, for which all tested models will be excluded by 

i 

the hypothesis of i being 'true'. Of course, such ID-value of M z > should be smaller (or 
at most equal), than the ILC discovery reach on model i. This procedure can finally 
be iterated, in turn, to all the different Z' models and the assessment of corresponding 
ID-reaches. This naive x 2 procedure can also be extended in a straightforward way to 
estimating exclusion ranges — and corresponding identification limits — in the cases where 
cos/3- and/or ctLR-dependent Z' models are considered in Eq. (jHJ). 

Examples of pairwise 'confusion' regions and corresponding contours, relevant to the 
Z' models chosen in Fig. [TJ are shown in Fig. [TDJ In this figure, the various steps of the 
procedure outlined above, as well as the final derivation of the ID-limits, can easily be 
followed. As an example of how to read this figure, consider the hypothesis that the rj 
model is 'true' (lower left panel), with M z > = 6 TeV. Then, if instead the ip or \ model 
should be true, the mass would have to exceed 4.2 or 6.3 TeV, respectively. 

Finally, Fig. [TT] shows the comparison of identification reaches or distinction bounds 
on the Z'-models considered in Fig. [TJ together with the corresponding bounds on M Z i 
obtained from the process pp — > l + l~ + X at the LHC with cm. energy 14 TeV and 
time-integrated luminosity 100 fb _1 . We assume, for the ILC, the same cm. energy, 
luminosty and beam polarization as in Fig. [IJ The figure speaks for itself, and in particular 
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Figure 10: Regions of confusion (yellow shaded areas) between a 'true' Z' model and the 
'tested' Z' models in the mass plane (M tTue , M tes ted) at 95% C.L. obtained from combina- 
tion of all polarized processes e + e~ — ► // at ^/s—0.5 TeV and £i nt =500 fb _1 . The dashed 
lines indicate the identification reach. 
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Figure 11: Comparison of the Z'-model distinction bounds on Mz> obtained from com- 
bined analysis of the unpolarized and polarized processes ([I]) at the ILC with -^i = 
0.5 TeV (1 TeV) and £ int = 500 fb _1 (1000 fb -1 ), compared to the results expected from 
Drell-Yan processes at the LHC at 95% C.L. [14]. Two options of polarization are con- 
sidered: unpolarized beams P = P + = and both beams are polarized, \P~\ = 0.8 and 
|P+| = 0.6. 

clearly exhibits the roles of the ILC parameters. In summary, one might be able to 
distinguish among the considered Z' models at 95% C.L. up to M z > — 3.1 TeV (4.0 TeV) 
for unpolarized (polarized) beams at the ILC (0.5 TeV) and 5.3 TeV (7.0 TeV) at the ILC 
(1 TeV), respectively. In particular, the figure explicitly manifests the substantial role of 
electron beam polarization in sharpening the identification reaches. Positron polarization 
can also give a considerable enhancement in this regard (if measurable with the same high 
accuracy as for electron polarization), although to a more limited extent in some cases. 

Clearly, our analysis is greatly simplified by the fact that the vector and axial vector 
couplings of the considered Z's are fixed theoretically. If we wanted to determine them in 
general, namely, with both masses and coupling constants a priori free variables, the \ 2 
analysis should be five- dimensional with, in addition, the limitation that for Mz> ^/s 
(contact-interaction regime), Mz> could not be simultaneously extracted. In principle, 
data at different collider energies could be utilized in this regard, for Z' masses not too 
far from [40]. 
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5 Concluding remarks 



We have explored in some detail how the Z' discovery reach at the ILC depends on the 
cm. energy, on the available polarization, as well as on the model actually realized in 
Nature. The lower part of this range, up to Mz> — 5 TeV, will also be covered by the 
LHC, but the identification reach at the LHC is only up to Mz> < 2.2 TeV. 

In this LHC discovery range, the cleaner ILC environment, together with the availabil- 
ity of beam polarization, allow for an identification of the particular Z' version realized. 
Actually, this ILC identification range extends considerably beyond the LHC discovery 
range. Specifically, the ILC with polarized beams at ^/s = 0.5 TeV and 1 TeV allows to 
identify all considered Z' bosons if Mz> ^ (6 — 7) x y^s. This represents a substantial 
extension of the the LHC reach. 
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